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A DESIGN PARAMETER SYNTHESIS DERIVED FROM A 


MATHEMATICAL ANALYSIS OF A HYPOTHETICAL 
LUNAR FLYING VEHICLE 

By William A. Jacqmein, David L. Hall, and Mark K, Craig 


SUMMARY 


The purpose of this study is to develop design information from a 
mathematical analysis of a manually controlled lunar flying vehicle 
found in NASA General Working Paper No. 10,079 (ref. l). The equations 
developed in that paper are discussed with respect to feasible vehicle 
configurations and the corresponding differences which arise from the 
‘use of thrust-vector as opposed to kinesthetic control. Experimental 
data concerning the human-pilot describing function are employed in the 
development of a stable region and a region of comfortable control 
expressed in terms of vehicle parameters. The equations representing 
vehicle parameters are investigated, and favorable lunar flying vehicle 
configurations indicated by this analysis are described. 


INTRODUCTION 


This paper is a study of the equations developed in the paper "A 
Mathematical Analysis of the Stability of a Lunar Flyirg Vehicle with 
Manual Controls.” In that paper, the vehicle equations of motion were 
developed using the Lagrange energy equation and were subsequently solved 
for the angular accelerations. The equations were then transformed to 
the state-variable matrix notation by use of the Jacobian. Vehicle con- 
trol was represented by the experimentally determined pilot describing 
function. The vehicle-control and -dynamics equations were combined; 
the closed- loop system characteristic equation was then derived and 
represented in Hurwitz matrix form yielding the principal minors. This 
study analyzes both the system characteristic equation and the system 
open-loop transfer function by means of various algebraic and geometric 
stability criteria to derive design-related vehicle parameter informa- 
tion. 


SYMBOLS 


vehicle paraneter term 

damping coefficient 
dissipation function 

derivative of ( ) with respect to q 

control force 

open-loop transfer function 
gravity 

Hurwitz-matrix principal minor 
moment of inertia 

pilot gain 
spring constant 

directed length from mass-1 center of gravity to gimbal point 

directed length from gimbal to mass-2 center of gravity 

directed length from gimbal to application of control force 

ratio of output to input 
mass 

general force term 
generalized coordinate 


radius of circle 



engine thrust 
kinetic energy function 

general lag-time constant 

general lead-time constant 

first-order neuromuscular system lag-time-constant approximation 
first-order neuromuscular system lag-time constant 

time 

potential energy function 

coordinate of mass n center of gravity 

body axis 

control-element transfer function 
pilot describing function 
coordinate of mass n center of gravity 
body axis 

constant term in vehicle parameter 

coordinate of mass-2 angular position 
gimbal angle 
matrix determinant 

damping ratio of the second-order component of the neuromuscular 
system 

coordinate of mass-1 angular position 



u 

X,X' constant term in vehicle parameter 

T pilot reaction time lag 
sin"^(|) 

w general frequency 

undamped natural frequency of the second-order component of the 
neuromuscular system 

-i -1 . partial derivative of ( ) with respect to q 

aq 

Subscripts : 

1 mass 1 

2 mass 2 

b body 

c control 

i matrix row designation 

J matrix column designation 

N neuromuscular system 

n integer 

p pilot 

Operators : 

( * ) first derivative of ( ) with respect to time 
(*') second derivative of ( ) with respect to time 

s Laplace transformation operator 
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SYSTEM PARAMETER ANALYSIS 


Discussion of Equations 

k 

The rotational stability characteristics of the proposed manually 
controlled lunar flying vehicle are described by the equation (ref, l) 


^ [l (>• - V22^l) " ^u] 

^ ^ s[^2a(’^L ■ a) ■ 2 *2 a]|® 

^ ^ - I)} ^ Vai = ° 

The formulation of this equation using vehicle parametef terms 
Ag^ and Agg represents a two-mass system manually controlled by 

thrust-vector manipulation about a gimbal point. The vehicle is con- 
sidered to be mass 1 and the engine, mass 2 (fig. l). If, on the other 
hand, the vehicle is considered to be kinesthetically controlled, the 
vehicle is designated as mass 2 and the pilot as mass 1, As a conse- 
quence, the pitch angle, which was 0 in the first case, becomes 3 
in the second case. The equation defining the rotational stability 
characteristics for this kinesthetically controlled vehicle can be found 
by examining the matrix describing vehicle rotational motion (ref. l) . 
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When the fourth line is rewritten ia terms of the gimhal angle 6, 
the representation of the pitch angle "becomes 

e - + A,^2® (3) 


Performing a Laplace transfozmatlon and rearranging tens yields 
the transfer function for the dynamics of a kinesthetically controlled 
vehi cle 
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The transfer function for the dynamics of a thrust-vector^-controlled 
vehicle, obtained through a similar procedure, is 



1 ^ 22 ^ 

6 ■ 2 
s 


(5) 


The equation defining the rotational-stability characteristics of a 
kinesthetically controlled vehicle is, therefore, that of a thrust- 
vector-controlled vehicle (eq. (l)), with the exception that is 

replaced by and is replaced by 


The A^j terms represent vehicle parameters extracted from the 

state-variable matrix describing vehicle rotational motion. These param- 
eters are (ref. l) 




-A23 = - i jk(M^ - Mg)[(M^ + Mg)!^ . M^M2L2(L^ . L^)] 

* * ^ 2 ^ 2 )! 

-^2U = - I “(«1 * «2)[(“i * “2)^2 
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( 9 ) 


(XO) 


The Hurwitz algorithm determines the criteria which the system 
characteristic equation (eq. (l)) must satisfy so that stability can be 
attained. Thu necessary and sufficient condition for equation (l) to be 
a Hurvitz polynomial (and, thus, indicative of system stability) is that 
the principal minors of the Hurwitz matrix all be positive. The prin- 
cipal minors of the Hurwitz matrix, as found in reference 1, ai’e 
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Another method of evaluating system stability employs the 
transfer fimction 

open-loop 

.... . 

' P B^TjS + i)(t„s + 1 ) (X 3 + 1 ) 

( 16 ) 
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stability in this instance is determined by replacing s in G(s) 
with J(i) and plotting the function G(jo)) for an increasing w. By 
the Nyq.uist criterion, the system is stable if the vector locus of 
G(J(o) drawn from w = 0 to w -► (in the sense of increasing la) 
encircles the point (-1,0) in a counterclockwise direction. If the sys 
tem is fouTid to be stable, its degree of stability can be approximated 
by constructing a circle with center on the negative real axis which is 
mutually tangent to the function G(jaj) and to a line from the origin 

at an angle where = sin" (l/M), The degree of stability can be 
evaluated using the equation for the radius of the circle 



(IT) 


An indication of the relative stability between two points can be 
gained by isolating the ratio of output to input M. An increasing M 
decreases the system stability. Acceptable stability requires that 
1 < M < 1.4 (ref. 2). 


Human-Pilot Describing Function 

The determination of the range of values assigned to the parameters 
of the human-pilot describing function is prerequisite to a numerical 
analysis of a man- vehicle system. The experimentally determined ranges 
of these parameters, as established by Skolnik (ref. 3), are the 
following : 


Gain , K 

P 

Lead time, Tj^ 

Lag time, T^ 

Reaction time, x 

Neuromuscular lag time, T^j 


0.6 to 250.0 
0.0 to 5.3 
0.0 to 25.0 
0.1 to 0.4 
: 0.1 to O.T 


It has been established by McRuer (ref. 4), however, that neuro- 
muscular lag time has little effect upon the pilot describing function. 
Consistent with this, McRuer demonstrates that a lag time of 0.1 yields 
results of good accuracy. 

Reaction time x is influenced by the complexity of the task, as 
well as by such tenuous parameters as fatigue, vibration, display 
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system, et cetera. Single-axis control experiments concerning human-pilot 
describing functions generally utilize a reaction time of approximately 

0.15. 


Experimental evidence indicates that the pilot is more satisfied 
with the control characteristics of a system when certain conditions 
have been met. At present, no precise data exist on exactly what 
constitutes this comfortable control system with respect to parameter 
definition. Various experiments, however, have yielded useful general 
observations. The most satisfactory system appears to be one in which 
the pilot acts as a simple gain and does not utilize the equalization 
constant, lead-lag compensation. Experiment 6 l 1 evidence indicates that 
the pilot prefers small values of lead, lag, and gain if this situation 
cannot be attained. Of these parameters, a small lead time appears to 
be the most decisive in the determination of what constitutes a com- 
fortable control system. According to these general observations, a 
comfortable control region can be approximated by limiting the range 
values of the pilot parameters to the following: 


Gain , K ... 

P 

Lead time, T_ 

Jj 

Lag time, 

Reaction time, t . . . . . 
Neuromuscular lag time, T^j 


0.6 to 5*0 
0.0 to 1.0 
0.0 to 5*0 
0.15 
0.10 


Development of Graphs 

Principal minors *bf the Hurwitz matrix (eqs. (ll) to (15)) can be 
used to develop maps of stable regions in the hold- 

ing the human-parameter values constant while selecting an arbitrary 
value of A^^ and subsequently solving for the value of A^g, which 

yields a minor value of zero. Trends established by all points gen- 
erated in this manner (fig. 2) determine the boundary between the 
positive and negative values of the minor and, thus, by the Hurwitz 
criterion, designate the periphery of stability in the 

Stability is attained in a region where each of the Hurwitz minors 
has a positive value associated with its 

(fig* 3)* The summation of each of the stable regions associated with 
a fixed set of human-parameter values over the entire spectrum of these 
values yields a stability map that is applicable to any given set of 
human parameters (fig. 4). A similar procedure is followed in 
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the development of a stability map for the special case of human- 
parameter values which are within the previously established comfortable 
region (fig. 5)* 


Redefinition and subsecLuent verification of these stability maps 
can be undertaken using the Nyquist criterion. Procedure dictates that 
the term s in the system open-loop transfer function G(s) (eq. (l6)) 
be replaced by jw and that G(jw) be plotted for an increasing o), 
holding all other terms constant (fig. 6). The previously discussed 
extension of the Nyquist procedure (eq. (iT)) provides a means of expres 
sing the ratio of output to input M as a function of the vehicle 

parameters (fig. ?)• In so doing, regions of stability defined by the 
ratio of output to input can be compared to those defined by the Hurwitz 
criterion (fig. 8). 


A qualitative indication of the interaction between the system- 
damping coefficient c and the system-spring constant k is presented 
in figure 9 for the specific lunar flying vehicle configuration depicted 
in figure 10. The interaction was developed by estimating vehicle char- 
acteristics corresponding to the stability boundary in the 

plane. Values of A^^ and A^g were then calculated using these char- 
acteristics and were subsequently set equal to the A. values of the 


The A.j equations 
Points were selected 


nearest point on the boundary of the stable region. 

(eqs. (6) to ( 10 )) were then solved for k and c. 
along other portions of the stability boundary and the process was re- 
peated. This procedure was followed for lunar flying vehicles (both 
kinesthetically and thrust-vector controlled) of various masses and iner- 
tias over the same range of A., values. 


SYSTEM CONFIGURATION ANALYSIS 


The state-variable-matrix equation describing vehicle rotational 
stability (eq. (2)) represents a gimbal-controlled two-mass lunar flying 
vehicle. Equation (2) imposes the constraint, however, that the centers 
of gravity of both masses be on the same side of the gimbal. Therefore, 
the applicability of this equation to cases in which the centers of 
gravity are on different sides of the gimbal point must be determined. 
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The state-VJiriable matrix is derived from application of Lagreuige's 
equation 


d f^^e\ 
^n ~ dt V j 


e ^ ^ 

3q^ 34^ 3q^ 


(18) 


The kinetic energy function is 


T 


1*2 *2 2 * 

[x^ 


•=i 

+ 2L 


«2 ^ 2a2 

0 + Lg 3 


cos 0 -■ sin 0^ 

+ 2L2$^X^ cos 3 - sin 3^ 


+ 2L^L203 cos(0 




^ 2 ^ 


(19) 


The dissipation function for rotational motion betwen t ¥0 bodies 


IS 


D =1 (e - B)' 


( 20 ) 


The potential energy function is expressed by 


V = M^gZj^ + M^g^Z^ + cos 0 + cos 3^ + ^ (0 - 3) 


( 21 ) 
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Referring to the kinetic energy term T (eq. (19))» the equations 
descrihing the velocity of mass 2 in relation to mass 1 are 

Xg » Xi + L^e cos 9-t cos 0 (22) 

Zg = - L^e sin 9 - Lg0 sin 3 (23) 


These are the critical relations in determining the applicability 
of the kinetic energy term to a configuration in which the center of 
gravity (c.g.) of one mass is above the gimbal point, and the c.g. of 
the other is below the gimbal point. In an arbitrary configuration 
(fig. 10), the vectors representing the motion of the gimbal point 
with respect to mass 1 and the motion of mass 2 with respect to the 
gimbal point are both negative in the X and Z directions . The second 
and third terms of Xg and Zg, therefore, should be negative. In 

figure 10, L^, Lg, 6, and 0 are also negative. Substitution into 

equations (22) and (23) yields 


Xg = X^ - L^e cos 0 - Lg0 cos 0 

(24) 

Zg = Z^ - L^0 sin 0 - Lg0 sin 0 

(25) 


Formulation of the relative-velocity equations from the configura- 
tion dynamics yields the same results. The relative-velocity equations 
(and hence, their conjugate kinetic energy term) are, therefore, in- 
dependent of two-mass c.g. distribution about the gimbal point. 

Investigation indicated that only the magnitude of the damping 
term D (eq. (20)) would be altered by the proposed configuration 
extension. Similarly, the potential energy relation V (eq. (2l)) 
need be considered only to the extent that the point of zero potential 
energy must be redefined to occur when and Lg are l80® apart. 
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The general force term Q- of the Lagrange energy equation 

n 

(eq. (18)) varies in form if the centers of gravity of the two masses 
are on different sides of the gimhal. Originally 


sin(6 - B) 

(26) 

= - F L 

2 c c 

(27) 

Q, = T sin 3 + F cos 6 
j c 

(28) 

= T cos 3 “ F^ sin 9 

(29) 


If, however, one center of gravity is above the gimbal point and 
the other below, the resulting terms are 


= - TL^ sin(0 - e) 


Q- = T sin 3 - F cos 9 
^3 c 

= T cos 3 + sin 0 


(30) 

(31) 

(32) 

(33) 


Examination of the equations describing vehicle rotational motion 
(eq. (2)), however, reveals that the F term affects only the con- 

stant added to the matrix and, hence, does not influence the system- 
characteristic equation. 

Thus, a term-by-term analysis of the Lagrange energy equation shows 
the system-characteristic equation to be independent of the distribution 
about the gimbal point of the centers of gravity of the two system masses. 
The characteristic equation can, therefore, be used Without restriction 
to indicate preferable vehicle configurations. 
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Stability maps (fig. generated by the Hurvitz criterion reveal 
that a necessary, but not sufficient, condition for the realization of 
vehicle stability is that parameter terms and have the 

same sign. Keeping this in mind, a useful generalization as to system 
preferability can be concluded from an examination of these A^^ terms 

(eq.s. (6) to (lO)). In simplified form they are 


A = 
21 

±ak ± X 

(3l*) 

II 

OJ 

CM 

< 

±ac 

(35) 

hi " 

ta'k ± X 

(36) 

h2 ” 

±a'c 

(37) 


The factors which influence the sign associated with an 
quantity are the sign and relative magnitude of and Lg. 


The 


system spring constant k is neglected, because a system in which the 
pilot will not be called upon to input a spring-constant value is 
preferred, A preferable vehicle configuration, therefore, requires 
that the signs of a and X be the same in order that the signs of 

A^^ and A^g be the same (a stable system). Figure 10 is an ex- 

ample of an unpreferable system. Assuming this configxaration to be 
kinesthetic ally controlled, the applicable vehicle parameter terms 
are A|^^ and A^^g. Mass 1 is the center of gravity of the man; 

mass 2 is the center of gravity of the vehicle, and the gimbal point 

is the man's ankles. Considering to be greater than Lg, which 

is normally the case, results in a system similar to current kines- 
thetic control vehicle simulators. When and Lg are substituted 

into equations (6) to (lO), the resulting relations are 



(38) 


* 1,2 = 


(39) 
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Because the signs of a and X are different, the pilot must 
input some spring-constant value to attain stability j consequently, the 
configuration is considered unpreferable . Figure 11 presents arrange- 
ments of centers of gravity and gimbal points which are considered 
preferable by the previously established criteria; some, of cotirse, 
are physically unrealistic. 


RESULTS 


Stability Maps 

Stability maps (fig. 2) describing the zero lines of the Hurwitz 
minors reveal, in accordance with theory, that minor is the deter- 
mining minor. If is positive, each of the other minors is also 

positive. The zero line of minor always passes through the origin 

in the plane. The stable region, therefore, does not transcend 

the abscissa and, thus, exists either entirely above or below the A.^ 

axis. Consistent with the positive minor stability criteria is the 
fact that the last term of the system-characteristic equation 

must be positive (K^ and must have the same sign). A negative 

K A^^ creates a root with positive real part in the solution, which 

is indicative of divergent and, thus, unstable motion. Since Ail and 

A ^2 must have the same sign, the vehicle parameters must place Ail 

and A ^2 values in the first quadrant for stability to be attained, 

if the pilot adopts a positive gain. If the pilot inputs a negative 

gain, on the other hand, stability requires that vehicle parameters 

A^, and A.^ be in the third quadrant. A vehicle whose A. param- 
il i2 ij 

eters place it in the second or fourth quadrant will be unstable. 
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Minor is .described by 


K 


H,. = 


H, 


(2 


' Vil 


2 i“3 2 


TjTn |a^2 + - g)J 


1^2 + T. 


2 

K ‘^A.. 

4- I. . P .. , } 


I *N 


]| 


^2 

hl[C^Z * -d 2* Wh - 1 ]} 


(14) 


The equation is unchanged if K^, '^il* ^i2 have the 

same sign. Because minor Hj^ is the critical term in stability deter- 
mination for both the first and third quadrants, the stable regions in 
each of these quadrants are of the same area (fig. 4 ). 

The comfortable region (fig. 5 ) is defined by the estimated region 
of comfortable control as derived from the previously developed restrained 
pilot parameters. The comfortable region is located near the origin and 
quite close to the approximate stability limit, rather than in the center 
of the stable area. An increase in the damping constant appears to be 
all that is necessary to vary the A^g parameter values from the sta- 
bility borderline tp the comfortable region. Although the magnitude of 
this increase varies from case to case, preliminary investigations indi- 
cate it to be large. 

A qualitative analysis of the interaction between the system-damping 
coefficient and the spring constant (fig*. 9) reveals that, as the system 
masses and inertias are increased, larger values of the damping coeffi- 
cient and spring constant are necessary to maintain system characteristics 
corresponding to a specific point in the A., plane. 

Ij 
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Stable-Area Dependence Upon Pilot-Describing Parameters 

The complexity of definition involved in vehicle-describing terms 
necessitates that only general observations be made concerning the 

dependence of stability maps upon pilot-describing parameters. The first 
of these parameters, pilot gain K inversely affects the stable region. 

Jr 

As predicted by. theory, an increase in the absolute value of the gain 
redaces the stable area. In a similar manner, an increase of the lag- 
time constant T^ also increases the area within' the stable border. The 

largest stable areas are approached as the lag-time constant approaches 
its upper limit of 25. An increase in the lead time tends to move 

the stable area toward the A^^ axis, thus reducing the range of A^g 
values. 


Nyq.uist Plots 

Redefinition of stability criteria in terms of Nyq.uist plots pro- 
vides a means of verifying results obtained by previous methods. Such 
a verification reveals the system to be conditionally stable; the gain 
can be increased or decreased to yield system instability. In accord 
with the previously established pilot-parameter range, the maximum gain 
allowing system stability is approximately 5 (fig. 6). 

The previously discussed extension of the Nyq.uist procedure 
(eq. (it)) provides a means of expressing the ratio of output to input 
M as a function of the A^^ vehicle parameters. Such an expression 

is presented in figure T for the pilot-preferred simple gain situation. 

A ratio of output to input in excess of 4 corresponds to the Hurwitz 
criterion A^j vehicle-parameter ste.bility borderline. Figure 8 is an 

extension and simplification of the correspondence between the stability 
limits determined by the ratio of output to input M and the stability 
limits detei'uined by the Hurwitz criterion. Since acceptable stability 
requires that 1 < M < 1.4 (ref. 2), the vehicle parameter-established 
stable regior contains only a small region in which acceptable stability 
is actually obtained. In the remaining area, the transient oscillatory 
response in the ratio of output to input is too extreme to be desirable. 
Thus, the majority of the area within the Hurwitz stability map has a 
very small degree of stability. 
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CONCLUSIONS 


Tentative conclusions, which are pertinent as parameters in the 
design of a lunar flying vehicle, can be drawn from this work. 

1. Vehicle-describing parameters and should have the 

same sign. For this reason, stable regions of eq^ual area exist in 
quadrants one and three of the 

2. The stable region is large to the extent that few restraints 
can be placed on vehicle characteristics without a knowledge of the 
system spring and damping constants. 

3. The degree of stability in much of the stable region is very 
small (fig. 8). 

k, A stable region exists where the pilot can act as a simple 

gain. 

5. The estimated region of comfortable control within the 
stability map is located close to both the origin and the stability- 
limit (fig. 5). 

6. The area of acceptable system response (fig. 8) indicates 
that a small value of the A^^ vehicle parameter is preferable. 

Because an increase in inertia of the controlled element decreases 
the value of A^^, large values of inertia are not necessariljr 

detrimental. Of primary importance, however, is the relationship 
between inertia and the other system parameters (mass, thrust, length). 

7. Stabilization of an unfavorable thrust-vector vehicle con- 
figuration requires smaller values of the system spring and damping 
constants than does stabilization of an unfavorable kinesthetic 
vehicle configuration. 

8. A kinesthetically controlled vehicle configuration allows 
less freedom in vehicle-design definition because the physical pilot 
parameters (mass, length, inertia) cannot be adjusted to meet stability 
criteria. 

9. Certain favorable vehicle configurations present themselves 

(fig. 11). 




REFERENCES 


1. Hall, David L.; A Mathematical Analysis of the Stability of a Lunar 

Flying Vehicle with Manual Controls. NASA General Working raper 
No. 10,079. 1968. 

2. D’Azzo, John J.; and Houpis, Constantine H. ; Feedback Control System 

Analysis and Synthesis. McGraw-Hill Book Co., Inc,, New York, I966. 

3. Skolnik, Alfred: Stability and Performance of Manned Control Systems. 

IEEE Transactions on Human Factors in ‘Electronics , Vol. HFE-7, 

No. 3, 1966, pp. 115-124. 

4. McRuer, Duane; Graham, Dunstan; Krendel, Ezra; and Reisener, Jr., 

William: Human Pilot Dynamics in Compensatory Systems. Air Force 

Flight Dynamics Laboratory, Technical Report No. AFFDLTR-65-15 * 

July 1965. 


BIBLIOGRAPHY 


Del Toro, Vincent; and Parker, Sydney R.: Principles of Control System 

Engineering. McGraw-Hill Book Co., Inc., New York, 19^0. 

Etkin, Bernard; Dynamics of Flight, Stability, and Control. John Wiley 
and Sons, Inc., New York, 1959. 

Jex, H. R. ; McDonnell, J. D. ; and Phatak, A. V. : A Critical Tracking 

Task for Manual Control Research. IEEE Transactions on Human Factors 
in Electronics, Vol. HFE-7, No. 4, I966, pp. 138-145. 

Summers, L. G. ; and Ziedman, K. : A Study of Control Methodology with 

Annotated Bibliography. NASA CR-125, 1984. 
















i 

1 


Figvire 3.- Stability interpretation of Hurwitz-minor zero lines 
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Figure U.- Stable-region definition in the A -A.^ plane. 




Approximate comfortable limit 


Gain, Kp - ± 1 .0 - i 
Lead, T. - 0 - 1.0 
Lag, T| - 0 - 5.0 
N.M. Lag, T^, - 0.1 
Reaction, T -0.15 


stability limit for comfortable human-parameter values in the A -A 


Figure 5 
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Figure 7.- 


Ratio M of output to input as a ftmction of vehicle parameters and 






Region of acceptable response 


Figure 8.- Comparison of regions as determined by the ratio 

M of output to input and the Hurwitz-minor zero line. 





















(t) Kinesthetic vehicle control 


